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simulation of the welding process, an experiment is designed 
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on a high strength steel cylinder, Fabrication of the HY-130 
cylinder, selection of welding parameters for the seam and 
girth welds, and the selection and arrangement of instrumenta- 
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element computer simulation of heat flow in weldments as well 
as techniques for axisymmetric finite element analysis of 
thermal stress and metal movement during welding are dis- 
cussed. Finite element representation of weldment geometry 
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in the form of photo-micrographs and electron scan micro- 
graphs. 
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CHAPTER I 


I.A INTRODUCTION: 

Welding has been used as the primary means of fabricating 
ocean structures for the past forty years. Design tools for 
the fabrication of welded structures have evolved primarily 
as the result of experimentation. These tools are usually 
expressed as empirical relationships among welding parameters, 
some of which include structure and weld joint geometry, 
physical properties of the structural materials and welding 
consummables, as well as methods and rates of energy trans- 
fer during the joining process. In most cases the use of 
these empirical relationships has been adequate to ensure the 
achievement of good quality weld joints in a production en- 


vironment. 


As design tools become more sophisticated, and the in- 
terrelationships among physical properties of welding and 
structural materials are better understood, design practices 
relating to allowable stresses in welded structures converge 
toward material limits. Development and refinement of tech- 
niques for detecting and classifying flaws in weldments is 
another design consideration which tends to increase confi- 
dence in welded structures subjected to more demanding re- 
quirements, allowing design engineers to increase allowable 


stresses or to reduce factors of safety or uncertainty. 
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An area of welding technology which requires further 
definition and refinement is the occurrence of residual 
stresses and structural distortion during and as the result 
of the welding process. Extreme non-uniform thermal cycles 
of weld metal and of structural metal in areas near the 
weldment produce unevenly distributed thermal stresses 
which cause localized plastic deformations, with resultant 
residual stresses and structural distortion. These residual 
stresses Can significantly contribute to structural failure, 
particularly when the failure mechanism is dependent upon 
low applied stress, such as brittle fracture or failure re- 
sulting from stress corrosion cracking. Distortions caused 
by welding can also cause structural failure by reducing a 


structure's ability to withstand buckling loads. 


The ability to accurately predict the magnitude and 
type of reSidual stresses and distortion resulting from the 
welding of ocean structures would greatly enhance a design 
engineer's ability to estimate total stresses which materials 
must withstand during loading. More accurate estimation of 
total stresses would allow for more effective use of materials 
and the increased predictability of structural failure would 
ultimately result in the production of more reliable ocean 


structures. 
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I.B BACKGROUND: 

Attempts to analytically predict or to simulate thermal 
stresses and metal movement which occur during the welding 
process have been the subject of considerable research effort 
during the past ten years. The complexity of the interrela- 
tionships among welding parameters which contribute to the 
formation of residual stresses in weldments does not allow 
the development of empirical relationships which are useful 
for more than one type of material, joint geometry, or welding 
process. Welding variables were too numerous and their in- 
terrelationships too complex for the development of a general 
technique for the analysis of transient thermal stresses, 
residual stresses, and distortion in weldments. Design 
engineers usually relied upon trial and error, or the use of 
design lanes derived from similar applications of proposed 
welding processes. 

Computer simulation of the welding process has offered 
reasonable opportunities to develop techniques which apply 
to wide ranges of materials, weld geometries, and processes. 
A description of the chronological and systematic development 
of techniques used for analysis of thermal stresses and metal 
movement during the welding process is contained in a paper 
presented by Masubuchi! to an international conference on 
Residual Stresses in Welded Construction and their Effects. 


This conference was held in London 15-17 November 1977, 
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hosted by the Welding Institute, Abington Hall, Abington, 
Cambridge, England. Systematic development of computer 
techniques for welding simulation began with a one-dimen- 
sional analysis by Tall in 1961. This research involved 
the simulation of thermal stresses occurring during the 
bead welding of a strip. Efforts at Massachusetts Institute 
of Technology by Masubuchi ” in 1968 expanded Tall's analysis 
into a one-dimensional FORTRAN program. This research was 
expanded into the development of a two-dimensional FORTRAN 
program in 1970 by ian Iwaki attempted to apply two- 
dimensional finite analysis to bead on plate welding. The 
most recent computer program to be used at MIT was developed 
by Muraki>. Muraki's program uses finite element analysis 
techniques to predict thermal stresses and metal movement 
based upon temperature distributions within the weldment 
which are calculated by heat flow programs also developed 


at MIT. 


fee SCOPE AND OBJECTIVES OF THIS STUDY: 

The purpose of this research is to design and carry 
out an experimental procedure to generate data which can 
be used to validate the existing heat flow program and the 
program which uses finite element analysis to predict the 
occurrence of thermal stresses and distortion during welding. 
The production of seam and girth welds in a high strength 


steel cylinder under very carefully controlled conditions 
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while measuring temperature, strain, and distortion as a 
function of time and location with respect to arc travel 
would serve as a valuable aid in evaluating the accuracy 

and potential usefulness of the heat flow program and of 
Muraki's axisymmetric finite element program for analysis 

of thermal stress and metal movement during welding. This 
research is part of a continuing effort at MIT to validate 
computer analysis of welding. A summary of past experimental 
studies of thermal stresses and metal movement during welding 
was presented by Masubuchi* and is reproduced and updated 


here as Table I-A. 
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CHAPTER ТТ 


MATERIAL CONSIDERATIONS 
II.A GENERAL: 

The Office of Naval Research is sponsoring an ongoing 
program to investigate thermal stresses and distortion in 
high strength steel (HY-130) structures. HY-130 steel is a 
high quality, low carbon, quenched and tempered steel with 
a yield strength of 130-150 ksi. This material was produced 
by the United States Steel Company, primarily for use by the 
Што Navy. 

ШЕВ COMPOSITION: 

Acceptable ranges of chemical composition for HY-130 
steel are shown in Table II-A. These ranges, reported by 
Shackleton^? are compositional limits of HY-130 steel as 
received for production. 

Mechanical properties of HY-130 steel in the "as received" 


quenched and tempered condition, as reported by Schroát?! are: 


yield Strength: 130-145 ksi 
Elongation in 2 inches: 15% 
Reduction of Area: 50% minimum transverse 


70% through thickness 
Impact requirements: 60 ft-lbf at 0° F 
A Rockwell hardness test conducted on the material to 


be used in this experimental study yielded a Rockwell "С" 


Scale Hardness of 32. 
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COMPOSITIONAL RANGES OF HY-130 STEEL Со 
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TABLE II-A 


. 60 
22/0) 
5 
. 40 
250 


05 


21710 


20220 


ОИ 


02:50 


02355 


51125 


07270, 


0.65 


O0 


.015 Мах 
Мах 
Мах 


.25 Max 





22 


ieee TEMPERATURE DEPENDENCE OF MATERIAL PROPERTIES: 


Subroutines within the axisymmetric finite element 
program developed by Muraki for analysis of thermal stresses 
and metal movement during welding require that material pro- 
perties be specified as a function of temperature. A search 
of the literature failed to produce detailed temperature 
relationships among the various material properties required 
by the Muraki program for HY-130 at elevated temperatures. 
Graphical presentation of these relationships reported by 
вснход+ 1 are reproduced in Figures II-l through II-10. А 
tabular summary of this information is included as Tables II-B 


and II-C, 


The Muraki program which is to be verified by experi- 
mental data assumes a linear relationship between various 
material properties and temperature. Assumed linear relation- 
ships which served as input data to the Muraki program are 


indicated by a dashed line on the appropriate figures. 





STRESS (ksi) 


150 


100 


50 


FIGURE IT-1 


25 


! - Compression 


- Tension 


0.004 0,008 9.092 
STRAIN (inches/inch) 


TYPICAL STRESS-STRAIN DIAGRAM FOR 
HY-130 AT ROOM TEMPERATURE 





YIELD STRESS (ksi) 


100 


50 





24 


0 600 1200 1800 2400 


TEMPERATURE (°F) 


FIGURE II-2 ESTIMATED EFFECT OF TEMPERATURE ON 0.2% 
OEESEL YIELD SGIRESS FOR HY-=-130 





) 


-3 


YOUNG'S MODULUS (ksi x 10 


25 


30 


20 


10 





0 600 1200 1800 8 2400 
TEMPERATURE (`F) 


FIGURE II-3 ESTIMATED EFFECT OF TEMPERATURE ON 
YOUNG's MODULUS FOR HY-130 





26 









-Compression 


- Linear Estimate 


150 
/ 


‘Tension 


100 


STRESS (ksi) 


50 


0 > 10 15 20 25 o 
TANGENT MODULUS (ksi x 10 


FIGURE ІІ-4  TANGENT MODULUS CURVES FOR HY-130 


3) 





TANGENT MODULUS (ksi х 1073) 


= 


30 


20 


10 


0 600 1200 1800 5 2400 
TEMPERATURE (TF) 


FIGURE II-5 ESTIMATED EFFECT OF TEMPERATURE ON 
TANGENT MODULUS FOR HY-130 
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SUMMARY OF MECHANICAL PROPERTIES 


ОРНОО ЕЕ 


ТЕМРЕКАТОКЕ 





а = Yield Stress 2% Offset (ksi) 
Е = Elastic Modulus (ksi х 107°) 
H = Tangent Modulus (ksi x 0 
V = Poisson's Ratio 


TABLE 11=В 
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SUMMARY OF PHYSICAL PROPERTIES 


FOR HY-1350 STEEL 





p = Density (1bm/in^) 


a,~= Coefficient of Thermal Expansion 
(microstrain/o,) 
C = Thermal Conductivity 5 
EIU ADEM CENE 
к = Specific Heat 
(BTU/lbm CF) 


TABLE II-C 
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CHAPTER III 
DESCRIPTION OF EXPERIMENTAL PROCEDURE 
ШЕР А FABRICATION OF THE WELD SPECIMEN: 


A circular cylinder of high strength steel was formed 
in order to test the MIT computer program for the analysis of 
heat flow in a weldment, with subsequent use of the tempera- 
Eure distributions from the heat flow program to validate the 
Muraki axisymmetric finite element program for analysis of 
thermal stress and metal movement during welding. The HY-130 
circular cylinder was formed to dimensions suggested in the 
June 1977 proposal on Study of Residual Stresses and Distor- 
Bon in Structural Weldments in High Strength Steels. This 
research proposal was made from Massachusetts Institute of 


Technology to the Office of Naval Research. 


Material for forming the cylinder was obtained from the 
Naval Ship Research and Development Center, Annapolis Labor- 
atory, Annapolis, Maryland. The material was received at 
Massachusetts Institute of Technology Welding Laboratory in 
the form of HY-130 plates 60" x 63" x 3/4". A request was 
made to Portsmouth Naval Shipyard in Portsmouth, New Hampshire, 
to roll the HY-130 plates into cylindrical shapes; Portsmouth 
Naval Shipyard being the only industrial activity in the 
Vicinity of MIT with experience in forming HY-130. Detailed 


drawings and specific requirements were submitted to Portsmouth 
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Naval Shipyard and an agreement was reached for the shipyard 
to form the plates into cylinders. However, due to the 
thickness (3/4") of the plate to be formed to a small outside 
diameter (18"), Portsmouth Naval Shipyard engineers suggested 
that it-would be necessary to form a cylinder by bending two 
plates to semi-circular cylindrical shapes and then welding 
the two halves together. In view of this requirement and in 
order to facilitate easier handling of the HY-130 plates, the 
two plates were flame cut to four plates 60" x 31%" x 3/4" 

at the MIT Welding Laboratory prior to transporting to 


Portsmouth, New Hampshire. 


Stages of cylinder fabrication are depicted in Figures 
111-1 and III-2. In order to prevent heat resulting from 
the flame cutting operation from affecting thermal properties 
of the steel close to the proposed weld, Tsai” suggested that 
flats required for bending be removed by flame cutting no 
closer than %" to the proposed leading edges of the semi- 
circular cylindrical mating surfaces. The flats were removed 
in this manner and remaining material to be removed in order 
to form a 60° Vee groove with 1/8" root, as illustrated in 
Figure III-2, was removed by mechanical grinding. Тһе semi- 
circular cylinders were then transported back to the MIT 


Welding Laboratory for experimental welding. 





STAGES OF CYLINDER FABRICATION 
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Stage 1 - HY-130 plate with dimensions show. 
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Stage 2 - Plate flame cut to dimensions shown 
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Stage 3 
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Plate halves rolled to semi- 
circular cylindrical shapes 
with flats approximately 1 3/)" 
remaining on leading edges. 


Flats renoved by flame cutting 
no closer than 1/l" to mating 
edges of semi-cylinders. 


l/l" excess material removed 
fron mating edges of semi- 

cylinders and 30° angle pro- 
duced by mechanical grinding. 
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ше В SELECTION OF WELDING CONSUMABLES: 


The type of welding suggested in the research proposal 
for girth welding of cylindrical shells was multi-pass gas 
metal arc welding (GMAW). This type of welding requires the 
use of shielding gas and filler wire, both of which can 
significantly affect weld quality. The most common type of 
shielding gas used in GMAW of low carbon steel is Argon with 
25% CO,. However, in the case of HY-130, experience has 


2 


shown that the use of Argon with 25% CO, leads to inadequate 


2 
fracture toughness properties in the weld deposit. Shackleton“ 
recommends rich Argon with less than 5% Oxygen or 10% СО, to 
avoid deterioration in the fracture toughness properties of 
the weld deposit caused by Oxygen pickup from the shielding 
gas. Shielding gas selected and purchased for this experi- 


ment was Argon with 1% O This gas was reported by 


27 
Radziminski~> as the gas most effective in shielding HY-130 


welds. 


Filler wire selected for the HY-130 cylindrical seam 
and girth welds was 1/16" diameter Linde 1405 wire. This 
wire was developed by Union Carbide specifically for HY-130 
weldments and was unanimously endorsed in the literature. Dur- 
ing periods when the filler wire was not in use, it was stored 
in a drying oven containing desiccative material. This step 
was viewed as necessary in order to prevent moisture or other 


impurities, or oxides, from collecting on the wire surface. 
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TIT.C SELECTION OF WELDING PARAMETERS: 


1. The variations in weldment geometry which could 
be considered for this experiment were somewhat restricted 
by the fact that cylinder size required that the weld be 
made only from the outside. Based upon several previous 
welding experiments reported by Doty^^ and the above men- 
tioned restriction, a single 60% Vee groove was selected as 
near optimum for producing the seam and girth welds of the 
cylinder. Weld joint geometry for this experiment is 


depicted below in Figure III-3. 


/ 


ve 60° Rn 
\included 
. angle 2 


А 









3/4" wall 
thickne 














| 8х" 


1/8" mating surface 


inside diameter 


FIGURE III-3 WELD JOINT GEOMETRY FOR PRODUCING 
SEAM WELDS ON HY-150%CYDLINDER 
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2. For the purposes of this study, heat input (Q) is 


defined in the conventional way as the product of arc voltage 
and welding current divided by arc travel speed, as expressed 
below in Equation 3-1: 


Nar 


HEAT INPUT = Q = р 





(Equation 3-1) 


There were some variations in welding literature regard- 
ing the upper level of heat input which would produce a good 
quality weld in 3/4" HY-130 steel. Shackleton*" recommended 
a maximum heat input of 40,000 Joules/inch. This upper limit 
was also recommended by Nunez” and Conner? Radziminski*” 
recommended a maximum of 45,000 Joules/inch. Although there 
was some variation in the suggested upper limits of heat 
input, there seemed to be general agreement that 40,000 Joules/ 
inch was a reasonable value for the production of good quality 
welds and this value was selected for production of the seam 
and girth welds on the HY-130 cylinder. Maximum and minimum 
allowable values of pre-heat and interpass temperature were 
also discussed in the above quoted references; general recom- 
mendations were more varied, but based on a rough average of 
recommended values, a common pre-heat and interpass tempera- 
ture of 200 °F was selected. 

3. Filler wire feed rate is automatically controlled 
by the GMAW welding machine. A wire feed rate of approxi- 


mately 87 inches/minute of the 1/16" Linde 140S filler wire 
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was measured during the production of a seam weld. The flow 


meter for the Argon/1%3 O, shielding gas was set at 40 et /hr. 


2 
A summary of welding parameters which remained unchanged dur- 
ing the course of this experiment is shown in Table III-A 
below. The only parameters which were varied between welding 
passes were the arc voltage, welding current, and in one 
instance the arc travel speed. Voltage, current, and arc 


travel speed were recorded for each pass and are included in 


the presentation of results. 


WELDING PARAMETERS 


JOINT GEOMETRY и ое 
Groove 


MAXIMUM HEAT INPUT 40,000 Joules/in 


FILLER WIRE 1/16" Diameter 
Linde 140S 


feed rate - 87 in/min 


SHIELDING GAS Argon/1% 0, 


flow rate - 40 fu isse 





TABLE III-A 
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III.D ARRANGEMENT OF APPARATUS: 


Welding equipment used for the entire experiment 

consisted of a Linde SVI-300 Power Supply with ST-2 Torch 

and associated arc travel and wire feed mechanisms. Seam 
welds were produced by placing the cylinder on an adjustable 
welding table and moving the welding torch along the groove 
of the seam weld. Prior to positioning the cylinder on the 
welding table, the two cylinder halves were tack welded at 
positions 6" from either end and at the center with 1" welds 
on the inside of each mating edge. Positions of the cylinder 
and welding torch for production of a seam weld are depicted 


in the photograph, Figure III-4, below. 





FIGURE III-4 POSITION OF CYLINDER AND WELDING TORCH FOR 
PRODUCTION OF SEAM WELD 
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ШЕРГЕ TESTING OF WELDING EQUIPMENT: 


In order to determine approximate temperature distri- 
butions during actual welding and to facilitate a thorough 
checkout of welding equipment prior to welding passes for 
data gathering, it was deemed advantageous to produce a seam 
weld on one side of the cylinder without benefit of instru- 
mentation to record temperature and strain. This allowed 
the wire feed mechanism, arc travel controls and indicators, 
gas flow regulators and meters, and the welding equipment 
electrical circuity and instrumentation to be tested and 
checked for accuracy under conditions expected during the 
actual test. This "test run" also provided an opportunity 
for operator training. A preliminary checkout of welding 
equipment had been conducted using "“bead-on-plate" welds on 


scrap pieces of HY-130. 


Approximate temperature distributions to be expected 
during the instrumented test were determined by the applica- 
tion of temperature indicating liquid (TEMPILAQ) in narrow 
bands transverse to the welding direction. TEMPILAQ liguids 


O O 


corresponding to temperature levels of 400 СЕ, 600 СЕ, 800 = 


Е, 
1000 Әр, and 1100 °F were applied to the chemically cleaned 
inner and outer surfaces of the cylinder.  TEMPILAQ was re- 
applied between passes and distances of temperature penetra- 


tion from the weld centerline on both the inside and outside 
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of the cylinder were recorded for each welding pass. Welding 
current, arc voltage, arc travel speed, and maximum distances 
from the weld centerline for the various TEMPILAQ levels were 
recorded for each welding pass and are shown in Table III-B. 
This tabulated data was used as a basis for determining strain 
gage and thermocouple locations for the instrumented seam and 
girth welds to be used in validating heat flow and thermal 


Frress computer programs. 
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WELDING PASS #* TEMPERATURE PENETRATION DISTANCE 


OUTS IDE INS DE 
LEFT RIGHT RIGHT 


-88 
40 
„35 


260A 
23V 
12 in/min 


350A 
22Y 
10 in/min 


3304 
217 
10 in/ min 


20 in/nin 


TEMPERATURE PENETRATION DISTANCES NOT RECORDED 


* Temperature penetration distances are expressed in inches, measured 
in the direction indicated from the weld centerline. 


ПЕВ ВЕБЕ = В 
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ІІІ.Е SELECTION AND ARRANGEMENT OF INSTRUMENTATION ТО 


MEASURE STRESS AND DISTORTION DURING WELDING: 


Stresses occurring during the thermal cycles associated 
with seam and girth welding of the HY-130 cylinder were de- 
termined by measuring strains as a function of time (i.e. - 
location of the strain measurement with respect to the time 
history of welding arc location) and as a function of distance 
from the weld centerline. Strains were measured by placing 
electrical-resistance wire gages on the outside surface of 
the cylinder. Bi-axial strain gages were placed with one 
axis aligned with the cylinder axis and welding direction 
(herein called the longitudinal direction) and the other axis 
aligned in a direction normal to the welding direction (herein 
referred to as the transverse direction. Characteristics of 
the strain gages used for this experiment are shown in 


Table III-C. 


In order to validate the MIT heat flow program, it was 
necessary to measure the heat distribution within the cylinder. 
Heat flow and temperature distribution within the HY-130 
cylinder was viewed as critical to validation of the Muraki 
program, inasmuch as strain, stress, and distortion calcula- 
tions are all based upon temperature distribution data which 
must be supplied as input to the program. Thermocouples used 


for measuring temperature distribution were Chromel/Alumel 
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STRAIN GAGE CHARACTERISTICS 


FAET-12D-12S6 (EL) 


Manufacturer: BLH Electronics, Inc. 





Waltham, Massachusetts 02154 


Resistance: 120 = .2 Ohms 


TABLE ILI=-C 


Strain indications were corrected for temperature 


by the following temperature error (apparent strain) equation: 
-2 2 -5m3 
pp EF 2. a T =a” a LOTA зи 10 "T 


- 2.40 х 1072т3 
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wires which were peened into small (0.040") indentations 
drilled into the cylinder outer surface. Thermocouple 
circuit continuity was assured by the use of an inductive 
spot welder on the thermocouple connections. After circuit 
continuity was assured, the thermocouples were protected by 


the use of sealing cement. 


Since longitudinal and transverse strains were to be 
measured with gages which required temperature correction, 
it was necessary to place thermocouples in positions which 
would give temperature-time histories consistent with strain 
gage locations. Arrangement of strain gages and thermo- 
couples on the HY-130 cylinder outer surface is shown in 
Figure III-5. Arranging the strain gages and thermocouples 
in this manner allowed them to be used for measuring strain 
and temperature distributions during the production of the 
second seam weld (the first seam weld was made without in- 
strumentation) and then to be used again after the cylinder 
had been cut into sections and mating edges prepared for a 
girth weld. Arranging the strain gages and thermocouples 
at a 45° angle with respect to the seam weld and at equi- 
distant locations with respect to each other allowed the 
strain gages and thermocouples to give strain-temperature 
distributions which could be adjusted in time so as to be 


Ecgncident. 
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It should be pointed out that the practice of locating 
strain gages and thermocouples so as to require time adjust- 
ment of their readings assumes quasi-steady-state heat flow 
(i.e. - the heat flow/temperature distribution with respect 
to the welding arc is the same at a position 40" along the 
length of the cylinder as at a position 41" along the cylin- 
der). This assumption is necessary in view of the fact that 
strain gages and thermocouples cannot be located at the same 


position. 


Arranging the strain gages and thermocouples at a 
45° angle with respect to the seam weld was viewed as ad- 
vantageous since this scheme places the strain gages and 
thermocouples which were subjected to the most severe thermal 
cycle during production of the seam weld in a position which 
1s subjected to the least severe thermal cycle during the 
girth weld. Failures occurring during the seam weld are 
most likely to occur at locations which are least critical 


to measurements to be taken during the girth weld. 


After determining and marking the desired strain gage 
locations, strain gages were cemented to the cylinder outer 
surface. Physical size and weight of the cylinder prevented 
the usual technigue of curing the strain gage bonding cement 
in laboratory ovens. In order to assure that the strain gage 
bonding would withstand expected temperatures, a 5 foot by 


> foot thermal insulating blanket (graciously manufactured 
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for this purpose by Bethlehem Steel Corporation Shipyard in 
Boston, Massachusetts) was placed around the cylinder and 
electrical heating elements placed inside. The temperature 

of the cylinder was slowly raised to 250 ОБ апа the strain 
gage bonding cement was cured at this temperature for a period 
of four hours, as directed by the strain gage manufacturer's 
instructions. The heating elements used for this purpose 


were later used to preheat the cylinder prior to welding. 


For production of the girth weld, it was necessary to 
cut the HY-130 cylinder into sections. Original plans called 
ШЕ спсбіпа the 60 inch cylinder into two 30 inch sections 
which would have mating edges prepared for a girth weld by 
machining a 609 single Vee groove. Due to physical limitations 
of laboratory facilities it was necessary to reduce the length 
of the two sections to be joined by a girth weld to 20 inches 
each. This required that the HY-130 cylinder be cut into 
three sections, each 20 inches in length. This decision was 
made prior to production of the instrumented seam weld, and 
the strain gages were placed so as to permit measurements 
during the seam weld as well as during production of the girth 


weld. 


After completion of the seam weld, the 60 inch long 
HY-130 cylinder was transported to the Army Materials and 
Mechanics Research Center in Watertown, Massachusetts, where 


it was sawed into three sections, each 20 inches in length. 





53 


A large power hack saw was used to cut the cylinder into 
sections. The saw was used on slow speed and the blade and 
cylinder were continuously cooled by an oil and water cool- 
ing mixture during cutting. This procedure was viewed as 
necessary in order to prevent heat generated by the cutting 
action from affecting thermal properties of the mating edges 


to be joined by girth welding. 


After the two saw cuts were made, producing the desired 
Бесе 20 inch long cylindrical sections, a ^ inch thick ring 
was sawed from the section which would not be used in the 
Girth weld; again paying particular attention to preventing 
heat buildup in the material being sawed. А diagram of the 


resulting sections appears below as Figure III-6. 





* Approximately 5 inch of cylinder length 
lost to width of saw cuts 


FIGURE III-6 
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A cross-section of the previously produced seam weld 
was removed from the ¥4 inch thick ring and later ground, 
polished, and etched for analysis of the metallurgical proper- 
ties of the seam weld. Details of this procedure and a dis- 
cussion of the results are contained in Chapter IV of this 


study. 


Mating edges of the two sections to be joined by girth 
welding were then prepared for welding by machining a 60° 
single Vee groove as shown below in Figure III-7. Heat 
buildup during machining was prevented by the use of a con- 


tinuously flowing oil and water cooling bath during machining. 


\ 
\ 60° : 
* included / 


d у апдје / 
ENT inch 
222 
1/8 inch 


, 


290 n. 20717. 
to cylinder end to cylinder end 





FIGURE Ill, 
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The two 20 inch cylindrical sections were transported 
back to the MIT welding laboratory where apparatus required 
for production of the girth weld was assembled. In order to 
control arc travel speed during production of the girth weld, 
it was necessary to design and manufacture a supporting struc- 
ture for the cylinder which would permit cylinder rotation at 
a controlled speed with respect to the fixed position of the 
GMAW torch. This was accomplished by producing a structure 
which was designed to support the two 20 inch cylinders from 
within. Details of this supporting structure are shown below 


in Figure III-8. 


+ ; быж 
| 1^ 3b 
| wm 
ji 
г. Ад а \ 


ту 
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The rectangular base of the supporting structure was 
bolted to the base of a tilting turntable and the 20 inch 
long cylinder halves were fixed to the support by force 
exerted on the interior walls of the cylinder by the threaded 
dowels. The turntable was then tilted so as to place the 
now 40 inch long cylinder (with 60% vee groove located in the 
center) in a horizontal position. The threaded dowels were 
adjusted in order to provide for constant clearance between 
the stationary GMAW torch and cylinder outer wall during 
cylinder rotation. Turntable controls were adjusted and 
calibrated in order to achieve desired arc travel speeds. 

A photograph of the apparatus in use during production of the 


girth weld appears below as Figure 111-9. 





FIGURE III-9 
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In order to measure radial distortion during production 
of the girth weld, two Bourns Linear Motion Potentiometers 
were placed at the inner surface of the cylinder. Character- 
istics of the linear motion potentiometers are shown below 


in Table III-D. 


ACCEPTANCE TEST RECORD 


FOR LINEAR MOTION POTENTIOMETERS 


TESTS PERFORMED SPEC LIMITS 
MECHANICAL TRAVEL ,........... 0.44 in max 
ЕТЕД ТОМ ЕОВСЕ .....:2........ ATOZ MIN 


NOTA ELEMENT RESISTANCE ..... 4750 to 5250 Ohms 


INSULATION RESISTANCE. ....... 50 Megohms min 
at 500 V DC 


NOISE OR CONTACT 


RESISTANCE VARIATION...... 100 ohms max 
END SETTINGS ................ 0.1 to 5.0% 
each end 
LINEARITY ......... we S * 1.08 VR max 


TABLE IIT 
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Displacement transducer circuitry and arrangement of 
linear motion potentiometers within the HY-130 cylinder are 
indicated in Figures III-10 and III-11 respectively. Prior 
to mounting the displacement transducers in position within 
the cylinder, a calibration test was conducted by deflecting 
the displacement transducers through measured distances and 
recording resultant voltages. A plot of this calibration 
test appears as Figure III-12. After completion of calibra- 
tion tests, the linear motion potentiometers were installed 
within the cylinder as indicated in Figure III-ll. Holes 
were drilled and threaded in the cylinder support mechanism 
and the linear motion potentiometer support rods were 
threaded and screwed into these holes at diametrically 
opposed locations within the cylinder. Potentiometers 
corresponding to Channels 1 and 2 of the recorder were 
located at positions 1/2 inch and 1 inch respectively in 


the cylinder axial direction from the weld centerline. 
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GOULD BRUSH 220 
TWO-CHANNEL RECORDER 





FIGURE III-10 DISPLACEMENT TRANSDUCER CIRCUITRY 


CHANNEL 1 DISPLACEMENT TRANSDUCER 
SENSING ELEMENT LOCATED 5 INCH 


de FROM WELD CENTERLINE.  INSULATED 
FROM CYLINDER INNER SURFACE 
BY CERAMIC BEAD. 


DISPLACEMENT TRANSDUCER SENSING E. 
MOUNTING RODS, SCREWED INTO 
CHEINDERZSUPBORT: 


CHANNEL 2 DISPLACEMENT TRANSDUCER 
SENSING ELEMENT LOCATED 1 INCH FROM 
WELD CENTERLINE. ALSO INSULATED. 


CYLINDER SUPPORT MECHANISM 





FIGURE III-ll DISPLACEMENT TRANSDUCER SENSING ELEMENT 
LOCATIONS 
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Strain and temperature were measured at the indicated 
locations during production of the instrumented seam weld; 
and strain, temperature, and distortion were measured during 
production of the girth weld. Equipment used to record the 
measured values included the following: 
Minneapolis-Honeywell H1508 24 Channel Visicorder 


Three Helland Universal Bridge Balancer Units, 
Model no. 82-6 (ser no's 1332,1307,1312) 


16 Channel Strain Gage Calibrator 


(one arm bridge adapter) 


Gould Brush 220 Two-Channel Recorder (distortion) 


General arrangement and operation of this equipment 
during production of the girth weld are indicated below in 


работе ТТТ-13. 





PEGURE-TIL=1S 
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EEG PRODUCTION OF THE INSTRUMENTED SEAM AND GIRTH WELDS: 
1. The Seam Weld: After the HY-130 cylinder had been 
suitably supported and aligned; instrumentation and other 
apparatus arranged, checked and calibrated, the instrumented 
seam weld was produced by making seven passes. Figure III-14 
below shows approximate dimensions and location of each pass. 
Values of weld parameters for each pass are tabulated in 


Table III-E. 


SEAM WELD PASSES 





FIGURE III-14 
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2. The Girth Weld: Production of the girth weld 
шоа Careful positioning of the rotating HY-130 cylinder 
with respect to the GMA welding torch. The girth weld required 
Six passes with approximate dimensions shown below in Figure 
III-15. Girth weld parameters are also included in Table 


METE. 


GIRTH WELD PASSES 





FIGURE III-15 
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SEAM WELD PASSES ес GIRTH WELD PASSES 
| WELD 
|, PASS NO. NO. 


uKÉR ITT ат јеси J weiss n ————— d —Á M —— P ——!———— — J—— '— ——————————————— n GOO 


I = Welding Current in Amperes 
= 


v = arc travel speed in inches/sec 


arc voltage in Volts 


TABLE III-E SEAM AND GIRTH WELD PASSES 
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СНАРТЕК ТУ 


METALLURGICAL CONSIDERATIONS 
IV.A PREPARATION OF SPECIMEN: 


When the 60 inch cylinder was sawed into sections, as 
previously discussed in Chapter III, a cross-section of the 
instrumented weld was cut from the 5 inch ring. The weld 
specimen was cut to approximately 3/4 inch by 3/4 inch and 
marked for identification as shown in Figure IV-l. Mounting 
for polishing was in thermo-setting bakelite resin, pressed 
into a 1% inch diameter by 3/4 inch long cylindrical shape 
under hydraulic pressure of 4200 psi at 250 Әр for 2 minutes 
using a Buehler press. 

The specimen was then rough polished by successive use 
ОҒ 60, 120, 180, 240, 320, and 400 grit 3M TRI-M-ITE Silicon 
Carbide paper. Fine polishing was accomplished by successive 
ШЕР ОР 1.0р, 0.3u, and 0.05u alumina polishing powder. Etching 
for photo-micrographs required about 15 seconds with 5% NITAL. 
The 500X photo-micrographs, Figures IV-2 through IV-5, were 
produced by the use of a LEITZ WETZLAR microscope with NPL 
50X/0.85 lens, using Polaroid type 52 POLAPAN film. Fine 
polishing was repeated and the specimen was etched using 1% 
NITAL prior to the production of 1200X and 6000X electron scan 
Micrographs. These electron scan micrographs appear as Figures 
IV-6 through IV-13. Locations of the photomicrographs and the 


electron scan micrographs are indicated on Figure IV-1. 
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© - Weld Metal 


(2) = Fusion Zone 
G) - Heat Affected Zone 
(7) - Parent Metal 


FIGURE IV-1 


DIAGRAM SHOWING LOCATIONS OF PHOTO-MICROGRAPHS 
AND ELECTRON SCAN MICROGRAPHS. 3/4" SAMPLE 
WAS CUT FROM THE SEAM WELD (SHOWING CROSS 
SECTION OF WELD - WITH NOTCH IN UPPER LEFT 
CORRESPONDING TO WELD CENTERLINE). AFTER 
POLISHING AND ETCHING, 500X PHOTO-MICROGRAPHS 
AND 1200X/6000X ELECTRON SCAN MICROGRAPHS WERE 
PRODUCED AT THE LOCATIONS INDICATED. 
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ELECTRON SCAN MICROGRAPH OF 
HY-130 WELD METAL. 1200X. 


13 NITAL. 





FIG IV-7 ELECTRON SCAN MICROGRAPH OF 
НУ-130 WELD METAL. 6000Х. 


1% NITAL. 
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FIG IV-8 ELECTRON SCAN MICROGRAPH OF 
HYy— EUSSONSZONS. 1200Х- 
13 NITAL. 
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FIG IV-10 ELECTRON SCAN MICROGRAPH OF 
HY 150 HAL: T2000. 13 NITAL: 
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FIG IV-11 ELECTRON SCAN MICROGRAPH OF 
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FIG IV-12 ELECTRON SCAN MICROGRAPH OF 
HY-130 PARENT METAL. 1200Х. 
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FIG IV-13 ELECTRON SCAN MICROGRAPH OF 
HY-130 PARENT METAL. 6000Х. 
1% NITAL. 
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CHAPTER V 


DESCRIPTION OF ANALYTICAL TOOLS AND PROCEDURES 


V.A INTRODUCTION: 

Analytical procedures used in this study to predict 
temperature distribution, strain, stress, and distortion 
occurring during production of a weldment were developed at 
Massachusetts Institute of Technology. Both the heat flow 
program and the program used to analyze thermal stress and 
metal movement use finite element techniques and are written 


in the FORTRAN-IV language. 


Discussion of these programs within this study will be 
limited to a brief synopsis of their theoretical background, 
presentation of assumptions and apvlicable equations upon 
which they are based, and a summary of welding parameters 
pertinent to production of the HY-130 weldment which were 
used as input to the programs. Detailed discussion of these 
programs, lists, card decks, and manuals describing their use 
are located within the Department of Ocean Engineering at 


Massachusetts Institute of Technology. 
V.B DESCRIPTION OF THE HEAT FLOW PROGRAM: 


1. Theoretical Background. The finite element computer 
program used in this study to generate temperature distribution 
in the HY-130 cylinder girth weld is Program B-4 in a report 


developed by Muraki*? and Masubuchi entitled "Computer Programs 
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Useful for the Analysis of Heat Flow in Weldments" at MIT in 


1974. This program develops temperature distribution based 
on heat conduction in quasi-stationary state from a moving 
point heat source in a plate of finite thickness. The 
solution assumes: 

а. The size of the heat source is small compared 
to that of the body under consideration. 

b. Heat is supplied instantaneously. 

c. Thermophysical properties do not depend upon 
ш- прегасиге. 

а. Heat is transferred in conduction only (other 
modes such as radiation and convection are not considered). 

е. The material of the body is uniform and iso- 
tropic in the thermophysical properties. 

The solution is based upon the following equation which 

expresses the difference between the temperature at time t 


(in seconds after the arc is initiated) anada the initial 


temperature: 
-y | -yr FUER -Vry 
T-T = — e 2K је 2K + E не жк | (Ес. 5-1) 
O TIC — | س د‎ = 
Е n=1 r т 
n n 
Where: с = х -~ vt 
e - < zs y? + 22 
ү „2 + y? + (2nh - =)” 
2 2 2 2 
ко = £ + y + (2nh + Z) 
а = heat intensity 
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2. Geometry of Weld Specimen and Description of 
Finite Element Representation of Temperature Dieser outed Om 
Program: 


Division of the HY-130 cylinder girth weld into elements 
FOr use in the heat conduction program is illustrated in 
Figure V-1. Node coordinates (x, y, and z values) for each 
of the 114 nodes are listed in Table V-A. The weldment is 
represented by finite elements which are arranged so as to 
Cause node locations for analytical determination of heat 
distribution to be coincident with locations of thermocouples 
which were actually placed on the cylinder outer surface in 


ВЕНЕ EO experimentally determine temperature distribution. 


Node coordinates listed in Table V-A are based upon ап 
origin (x=y=z=0) which is coincident with the starting point 
of the heat source (welding arc), with subsequent times of 
interest corresponding to location of the heat source (in 
nodal coordinates) at the referenced time. This analysis 
assumes that the heat source moves on the plane z=0 along 
the x-axis. The x direction (direction of welding arc 
travel) would appear on Figure V-1 as normal to the y-z 
plane, positive direction into the paper. Іп this figure 
the y-axis is along the cylinder length and the z-axis is 


in the radial direction. 
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In order to determine temperature distribution at 
weldment nodes for the first pass of the HY-130 cylinder 
girth weld, the following parameters were used as input to 


the heat flow program: 


Velocity of welding arc (v) = 0.2000 in/sec 


Heat intensity (q) = 6084 Joules/ ас на 
Where: q = L x VOX Nare Ое o2) 
h 
I = welding current 
V = arc voltage 
h = plate thickness 
N gr welding arc efficiency 
Density (p) = 0.282 lbm/, 3 


Specific heat (x) = 124.4664 Joules/ + tbm Ор 


Thermal conductivity (с) - 0.5250 Јоч1е/ ec in Ор 


It should be noted that this stage of the computer 
analysis does not provide for temperature dependence of 
material properties. Values used in this case to determine 
temperature distribution are those valid at 400 Ор, It is 


also significant to note that for this analysis, an arc 


efficiency CD. of 0.65 was assumed. 
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No modifications were required to the analytical or 
mathematical techniques required for generating temperature 
distributions to be used by the Muraki strain, stress, and 
distortion calculating program. However, since large amounts 
of numerical data were generated by the heat flow program for 
input to the Muraki program, it was necessary to modify the 
heat program output. Since weldment geometry and finite 
element representation of the nodes and elements were the 
same (except for coordinates) for both programs, the heat flow 
program was modified so as to provide its output in the form 
of punched cards suitable for use as input to the Muraki 
program. This reduced the probability of introducing numeri- 
cal error between the two programs as well as improving the 


ease of program interaction. 
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WC” DESCRIPTION OF THE PROGRAM TO CALCULATE STRAIN, STRESS, 


AND DISTORTION: 


1. Theoretical Background: The program used in this 


study to analytically generate strain, stress, and distortion 
distributions within the НУ-130 cylinder girth weld was de- 
veloped at MIT by Toyohiki Muraki for thermo elasto-plastic 
analyses of axisymmetric structural components. The Muraki 
program is written in FORTRAN IV language and uses finite 
element techniques in an incremental fashion to calculate 


strains, stresses, and distortion. 


The Muraki program has great flexibility in terms of 
allowing temperature dependency of material properties to be 
considered as well as considering the effect of melting and 
deposition of metal which occurs during welding. Weldment 
restraint and concentrated loads are also considered in the 
analytical solution. In the "Manual on Axisymmetric Finite 
Element Program for Analysis of Thermal Stress and Metal 
Movement During Welding", Muraki^? presents a brief synopsis 
of the theoretical background upon which this analysis is 


based. This synopsis is reproduced here as follows: 


Deme role of Virtual Work 


The principle of virtual work for the quasi-static 


problem is, in general, 


ил £68) olav - јлл (бај СрЈам – јла (бај (F)ds - 0 (1) 
J 





82 


where superscript T stands for the transpose of a vector or 
a matrix. te. (0), and (а) are vectors of the time rates of 


strain, stress and displacement, respectively. 


° ° 


P and F denote the time rates of body force over the 
volume, V, and surface force on the boundary surface, 51, 
respectively. 


The boundary condition on the specified displacement on 


the surface, S is given as 


2: 
Шигг- (0 ол 5. (2) 
It is noted that the problem here is formulated in the 


small displacement theory based on the rectangular cartesian 


coordinate system. 


Constitutive Equations 


One of the important features in the analysis of thermal 
Stress and distortion during welding is the temperature de- 
pendency of material properties and the yield criterion. 
Therefore the constitutive equations are provided here in a 


general form. 
The rate of strain for the elasto-plastic state is given 
by: 


р = (Е у + {Е} (3) 
where superscripts e and p refer to elastic and plastic com- 


ponents of the strain, respectively. 





22 
The rate of plastic strain is assumed as 
(c P) - Ма) (4) 


where A is a proportional constant and g is derivatives of 


EE vlıeld function, f, with respect to the stress, o. 


12 
ШЫ - 1,2,3); 
Ir: Е 
98. (5) 
Using Hooke's law and Equation (4), we have 
• = o ж 
(ë) = =|De|~* {6} + Ag} (6) 


where E|De| is the matrix form of Hooke's law. 


Thermal strain and the temperature dependency of 


Material properties are added to Equation (6) as follows: 


: — b = © 2 == 
1€) 9 Ере| 1 160- Е, |ое| {о} | 
А (7 
n СА + Мај 


where HE is the thermal strain, and 5 is the Kronecker 


symbol. 
Equation (7) is rewritten for the rate of stress as, 


(о)- Е|ре| (ë)+ #(o) - е E|De| (6,.) 


- ДЕ пе| Та) (8) 


ж = 
Note: | | : stands for the inverse of the matrix 
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If the von Mises' yield criterion is adopted, £ in 
Equation (5) is given by 
Е=<6 - с (e, mq (9) 
13 
where C is related to the strain hardening and given as a 
function of plastic strain, є (В) апа temperature, T 
0 15 referred to as the equivalent stress and given 
by 
s цар ча (10) 


Te is noted that (c^) is the vector of deviator stresses 


which is different from іс). 


When a state of stress in a portion of a body is plastic 


and then changes to a new state of stress which is still 


plastic, the following relation is obtained from Equation 


(9): 


a fa 
aj jhh 
al 
+ 
rn 
ده‎ 
Fh 
= 
H 
r 
ro 
to 


ЗЕ ST л) (11) 


Equation (11) with Equations (4), (9), and (10) leads 


to 
^ Ue = > of s 
{o°} {o} = -2/30(- AH” + E T) (12) 
where 
2 ОЁ 
Н” = - —— 
and ZU! = y 2/3136 P!) }T{3e P) } 


rom Equations (8) and (12), A is obtained as follows: 
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AAA eto ne E) 


2/3 5 4" + Elo”) *|De| 19) 


E E (0935 (0) - аа леа 
Е 13 
ШЕ С 25 T) (13) 


Substituting A of Equation into (8) we have a general 


form of the stress-strain relationship as 


(9) = (E|De| - 1/S E|De|(g) Etc^)" |De|) x 
cs £^ (84.1) 
+ E/E (|1] - 1/5, E|Del (gH (o^)7) tc) 
ЗЕ 


= 1/5. Е|Ве| (а) 2/30 - Т 


(14) 


where 
S. - 2/3 v H^ + Elo") |ре| (а) 
and |I| is called the unit matrix. 
Equation (14) covers the stress-strain relationship in 
the elastic as well as plastic state. This can be easily 


confirmed by eliminating terms related to the yield function 


from Equation (14). 
Yield Criterion 


As described in Equation (9) the yield function used 
here is that proposed by von Mises. 
The yield function is to separate the region of elastic 


deformation from that of plastic deformation. Since the 
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deformation is closely related to the loading process, the 
yield function is also called a loading surface which is 


defined in stress space. 


If a part of a material in question is in the region of 
plastic deformation, which means that the state of stress of 
the part is located on the loading surface, one of the 


following states will be realized at the next instance: 


мола < о unloading) 
Ё = 0, f“ = 0 (neutral) 
ЕАО, Е > 0 (loading) 
where 
_ of ЭЕ 
ел Dij + эт т 


13 
If a part of the material is in the region of elastic 
deformation, one of the following states will be realized 


at the next instance: 


SO (elastic) 
(46) 
к= 0б (plastic) 


TES worth looking into Equation (15) more closely. 
Using Equations (10) and (12) with Equation (15), the 
following relationship is obtained: 

f = AH (17) 


Equation (17) indicates that when non-strain hardening 


is assumed of a material, Equation (15) cannot be used to 
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judge whether or not a part of the material recovers from 
the plastic to the elastic state of stress. To avoid this 


Ami culty, Equation (17) will be written in the form of: 


lim £- и 
Нс+0 22 : (18) 


Therefore Equation (15) is written as, 


f= 0, A 


A 


0 (unloading) 


ШИ УРО ЛД ОШ пеп тл) ШЕ) 
l= O, 2010-0106) 


to cover non-strain hardening materials. 


Finite Element Approach 


Various types of axisymmetric finite elements can be 
derived based on Equation (1). 


In the computer program described in the Muraki^? 


Manual, a four node-isoparametric element is used.  Follow- 
ing a standard procedure to derive an element stiffness 
matrix and element forces, the coordinates and the displace- 


ments for the elements are assumed by 


=1 11 
(20) 
4 
4 
ныг. (Ab 
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4 
DN Uh V (21) 
1 


where h ћ ћ and û, are interpolation runctions, х, 


Ш 2:3, 4 
апа are nodal coordinates, and u; and V; are nodal 


displacements. 


A key of the above derivation is use of the same in- 
terpolation functions for both the element coordinates and 
the element displacements. Details of the interpolation 


Eimetions are 


(22) 
йй = (1 + r)(l +s) 


Һ = 5(1 — r)(1 + s) 


where xr and s are referred to as natural coordinates and vary 


from -1 to +1. 


The element strains in the axisymmetric condition are 


€ да 
P 3x 
ду 
е — 
Y y 
= (23) 
u 
с. = 
23 
Y шэнэ 
ху 37y + 
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Because u and v in the strain are a function of x and 
y, the following relation is first required between the x-y 


and r-s systems 


9_ 9х ¥ д. 

or or r ох 
= (24) 

9 9х ду 9. 

ЕС | Js 55 | Jy 


Derivatives of the displacements with regards to r and 


Б ате 

ju x oh. 

eee sr u; 
au 

E la 

95 TELE 1 225) 

ШҮ: 

or 1=19Ё 1 

” х В. 

= 7 š Уі 
mu 


Using Equations (24) and (25), the strains in Equation 


(23) аге given in terms of nodal displacements, u u u 


DOT 
and u) as follows: 


(£) - |B| fq) (26) 





90 


where 


217 ze , 220) Est Y } 


оту : : : ! : 
(ај = fu 22027720 уд} 


and {B} is the matrix which gives a relationship between 
the strains and the nodal displacements. It is noted that 


ШИЕ 15 а function of r and s. 


The stress-strain relationship of Equation (14) is 


rewritten in a simple form as: 
(о) = |D| t£) 4 (6?) (27) 
where 


° ДЕ = e e “ 
{o} = {or 2111 20) 


{D} is the stress-strain relationship matrix and {5°} is the 
stress vector due to temperature change, and the temperature 


dependency of material properties and the yield function. 


The element stiffness matrix {k} is obtained using the 
principle of virtual work of Equation (1) and Castigliano's 


theorem as, 


{к} = S//1B)*1D)1B)y dodxdy (28) 


27/ (в) 1 (п! (в) у | 1 аказ 
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where 
ox эу 
Or dr 
|| = 
ax Эу 
95 95 


Introducing a numerical integration, we have 


T 
[k) = 212 Св); 510108}, 4191 (29) 


373-0; - 
i,j 1411 17 


where subscripts i and j indicate that Yas: {B}, and |J| are 
evaluated at points i and j. ЛЕ is weight factor at points 
i and j in the numerical integration. 

The program in the Muraki manual uses the Gauss-Legendre 


numerical integration and Е in that scheme can be found in 


text books on numerical integration. 


The element force vector Ta corresponding to the second 
term in the right side of Equation (27) is obtained in a 


similar manner. 


(#9) = 55148) (09 )удвахау 


ee (30) 
anss{B} {o }y|aldrds 


Again introducing the numerical integration we have 


6 E 
{Е} = 272 о ?Yis ld] (31) 


1491: 
1.3 го 


In the current version of the computer program described 


in the Muraki manual, only the static load considered is the 
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is the concentrated load at nodal points. Therefore there 


им other element force. 


The stiffness matrix and force vector for the whole 
body (i.e. the whole system) are obtained by assembling the 


element stiffness matrices and force vectors. 


The equilibrium equations for the whole system are 
obtained as 


{kJ}, {QO} + {f°}, = {Р} (32) 


т 1 J J 


where i and j stand for the i-th element and the j-th node 


respectively. 


{Q} is the vector of the time rate of nodal displacement. 


{F} is the vector of the time rate of concentrated load. 


Equation (32) is further written in the incremental 


form as follows: 


6 en 
| o xb Е J; “ийг (32) 


where 
a = oe 
20 
ue 


At is the time increment. 


ТЕ is noted that the time increment is included in 


ма лоп (33) implicitly. 
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Once the nodal displacements for the system, {AQ}, are 
obtained, the element strains are calculated by Equation (26) 


and then the corresponding stresses by Equation (27). 


When there is no further increment of concentrated and/ 


or thermal loads, the analysis is complete. 


When there is another increment of the loads, new 
equilibrium equations are calculated using the nodal dis- 
placements and stresses obtained in the preceding calculation. 
This procedure of incremental analysis continues until the 


Ши increment is solved. 
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2. Geometry of Weld Specimen and Finite Element 
Representation of Weldment used by Muraki Program: 


Finite element representation of the weldment is 
depicted in Figure V-2. A comparison of Figure V-2 with 
the finite element representation of the weldment, Figure V-1, 
used in the heat flow program will show that an axis trans- 
formation as well as a shift of the origin were required in 
order to use heat flow program output as input for the Muraki 
Program. 

For the Muraki program, the y-axis corresponds to the 
radial direction and the x-axis is along the cylinder length. 
The z-axis corresponds to the direction of welding arc travel, 
normal to the x-y plane as shown in Figure V-2, positive 
direction out of the paper. Тһе origin (x=y=0) for the Muraki 
program analysis corresponds to a point at the center of the 
cylinder end. z=0 corresponds to the time of arc initiation. 

The Muraki program also required the specification of 
nodal coordinates (x,y values). Element connectivity was input 
to the Muraki program by identification of the nodes, in counter- 
clockwise fashion beginning with the node closest to the heat 
source, associated with each element. Nodal coordinates input 
to the Muraki program are listed in Table V-B and element con- 
nectivity is shown in Table V-C. 

Thermophysical material properties required by the 


Muraki program were taken from Figures II-1 through II-10 
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in Chapter II of this study. Welding parameters required as 
input data for the first pass of the girth weld were the same 
as those used by the heat flow program. These parameters 


were discussed and values presented earlier in this chapter. 
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TABLE V-C (Cont'd) 


TO 
СНАРТЕК УТ 


Comparison of temperatures and values of strain, stress, 
ВОО 5ЕОГЕЗОП during production of the HY-130 cylinder girth 
weld with those predicted by use of the heat flow and Muraki 
programs follows in the form of tabulated data and in graph- 
ical presentations. During production of the girth weld, arc 
travel began at the seam weld opposite the cylinder seam about 
which thermocouples and strain gages were placed. Times at 
which data comparisons are presented were selected as repre- 
senting a thermal cycle during which temperature, strains, 
stress, and distortion underwent greatest changes. 

Temperatures determined experimentally at the eight 
thermocouple locations for the times of interest are shown 
in Table VI-A. Distances from weld centerline are indicated 
beside the column headings designated T1, T2, T3, T4, T4', 
T3', T2', апа ТІ". These temperatures are presented graph- 
ically in Figure УІ-1. Temperature distribution predicted by 
use of the heat flow program is presented by a computer gener- 
ated graph, Figure VI-2. 

Qualitative comparison of the analytical and measured 
temperature distributions shows remarkable coincidence. The 
slightly higher overall values of the analytically generated 
values may be explained by the assumed welding arc efficiency 


EB 2) эг 0.65. 
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For the purpose of presenting comparative strain, stress, 


and distortion data, the following terms are defined: 


SA) 


ВИЕ A 


се я 


А 


ВЕ) Хх 


BST ER) A 


EXPERIMENTALLY DETERMINED STRAIN 
(MICROSTRAIN) AT A POINT *.* INCHES 
FROM THE WELD CENTERLINE. 


ANALYTICALLY PREDICTED STRAIN 
(MICROSTRAIN) AT A POINT 7.3 INCHES 
EROM THE WELD CENTERLINE. 


EXPERIMENTALLY DETERMINED STRESS (1bf/, 2) 
AT A POINT *,* INCHES FROM THE WELD 
CENTERLINE. 


ANALYTICALLY PREDICTED STRESS (1b£/, 2) 
AT A POINT PENCHES FROM THE WELD 
CENTERLINE, 


EXPERIMENTALLY MEASURED RADIAL DISTORTION 
IN INCHES AT 2A POINT M: INCHES FROM THE 
WELD CENTERLINE. 


ANALYTICALLY PREDICTED RADIAL DISTORTION 
НЕ Е ПРО SE INCHES FROM THE 
WELD GElNTERLINE 


Positive values of strain and stress correspond to the 


tensile condition and negative values correspond to compressive 


SEFalns and stresses. 


Positive values of radial distortion 


correspond to displacement in an outward direction from the 


cylinder center. Tables and graphs on the following pages 


use the terms as defined above to compare data. Strain and 


Stress directions correspond to directions shown in Figure V-2. 
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Tables VI-B through VI-E present in tabular form 


comparisons between strains measured during production of 

the HY-130 cylinder girth weld and those predicted by use 

of the Muraki program. Tables VI-F through VI-I present 

the analytically predicted stresses. Experimentally deter- 
mined stresses were not tabulated since experimentally 
determined strains could not be transformed into stresses. 

There are three (ignoring creep) components of strain; elastic, 
plastic, and thermal. Elastic and thermal constitutive rela- 
tions are simple and the stresses resulting therefrom could 

be calculated if these strains could be separated into their 
components. However, the constitutive relationships among 
elastic, plastic, and thermal strains are complex, made further 
so by the presence of a strain hardening material and by the 
dependence of all material properties on temperature. For 

this reason, experimentally determined strains could not be 
transformed into stresses. Values of radial distortion measured 
during production of the girth weld and predicted by the Muraki 
program are shown in Table VI-J. Comparative graphs of measured 


and predicted data are shown in Figures VI-3 through VI-9. 
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LONGITUDINAL STRAIN TRANSVERSE STRAIN 
(X-DIRECTION) (Z-DIRECTION) 


STR(1.0)X_| STR(1.0)A STR(L.O)A 


55. 0 420 
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I5. 0 


7.5 2998 
TE = 
р 1539 Сен 





160.0 1930 | am | 680 | 
170.0 L 2200 : 714 
| 2000 | me | 290 | -4 726 


TABLE VI-B EXPERIMENTAL AND ANALYTICAL STRAIN 1.0 IN FROM 
WELD CENTERLINE 
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LONGITUDINAL STRAIN TRANSVERSE STRAIN 
(X DIRECTION) (Z-DIRECTION 


Шин Шї ЗБ) А | STRO EIX 


musasasa 
O olam 
O „. г ole 
LAE a 

EN 

u. 




















o. aml e | a 


55.0 -400 5.2 / 513 


| 160.0 — - минь ос а m= — нь - nm | - = er = e ————À 


A 


150.0  '| et 651 427 530 





TABLE VI-C EXPERIMENTAL AND ANALYTICAL STRAIN 1.5 IN FROM 
WELD CENTERLINE 
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LONGITUDINAL STRAIN TRANSVERSE STRAIN 
(X-DIRECTION) (Z-DIRECTION) 


EE. | STR(2.0)X ОӨТБ( 20012 STR(2.0)X STR(2.0)A 


| ЕСІН ШЕГЕН ШЕГЕ БЕТІН ШЕГЕ 
ата | ње | о | = Í ol 
ues | ње БЕРЕГЕН a 2 
Dumas | m T S T a L me 
E 


Ds [oe [ae 
ere 
ЕЕ sS ual wm] mm 
БЕТТЕН БЕ” ЕН БЕТТЕН БЕГЕН ШЕКТЕН 
ва | ae гэн Гэв ыг 
Due [um a ыг 

EN | es oe L me 
ues [oum | |І ee | 


150.0 | -384 | -216 606 378 | 


























TABLE VI-D EXPERIMENTAL AND ANALYTICAL STRAIN 2.0 IN 
FROM WELD CENTERLINE 
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LONGITUDINAL STRAIN TRANSVERSE STRAIN 
(X-DIRECTION) (Z-DIRECTION) 


SI A ОТВЕЛА ОДЕ ЕЈ SUR (2.5) A 


Еј 
EE | — |e | ms | 
Г Гош | - 
ЕГІХЕН БЕРЕН БЕРЕГЕН ШЕТЕН ШЕНГЕН 
Dues | ue (| — |o њо 
ume | ue | ow | зе [ш 
Dues | | | 39 


40 
0.0 2125 -114 340 | 205 


, 145.0 EMILE -123 355 244 


| 150.0 -249 Ae) 457 | 276 





200.0 = 176 203 365 463 


TABLE VI-E EXPERIMENTAL AND ANALYTICAL STRAIN 2.5 IN 
FROM WELD CENTERLINE 
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TABLE VI-F PREDICTED X AND Z DIRECTION STRESSES 1.0 IN 
FROM WELD CENTERLINE 
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LONGITUDINAL STRESS TRANSVERSE STRESS 
(X-DIRECTION) (Z-DIRECTION) 
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TABLE VI-G PREDICTED X AND Z DIRECTION STRESSES 1.5 IN 
FROM WELD CENTERLINE 
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LONGITUDINAL STRESS TRANSVERSE STRESS 
(X-DIRECTION) (Z-DIRECTION) 
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TABLE VI-H PREDICTED X AND Z DIRECTION STRESSES 2,0 ІМ 
FROM WELD CENTERLINE 
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LONGITUDINAL STRESS TRANSVERSE STRESS 
(X-DIRECTION) (Z-DIRECTION) 


STS (2.5)X | STS(2.5)A STS(2.5)X 
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TABLE VI-I PREDICTED X AND Z DIRECTION STRESSES 2.5 IN 
FROM WELD CENTERLINE 
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RADIAL DISTORTION RADIAL DISTORTION 
(Y-DIRECTION) (Y-DIRECTION) 


SON DSTO SEA ОБОДА став ОА 


—— о 


25.0 0.0000 
0.0000 0.0002 0.0000 


137.5 0.002 00032 0.002 0.0039 























0.0000 0.0000 


0.0000 0.0000 


0.0000 0.0000 


0.0000 0.0000 


0.0000 


140.0 0.003 0.0035 0.005 0.0041 


*DEFLECTION TOO SMALL TO QUANTIFY 


TABLE VI-J EXPERIMENTAL AND ANALYTICAL DISTORTION 0.5 IN AND 
1.0 IN FROM WELD CENTERLINE 
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CHAPTER VII 


CONCLUSIONS 


This study of residual stresses and distortion in high 
strength steel weldments has shown that existing computer 
programs require further refinement in order to be used as a 


day-to-day tool of the welding engineer. 


The occurrence of residual stresses and distortion during 
production of a weldment is a complex phenomena. Past efforts 
to simulate the welding process involved clever simplifying 
assumptions followed by attempts to carry out experiments which 
supported these assumptions. The new approach to welding 
simulation employs the enormous power of computer techniques 
to analyze what actually takes place during the welding process. 
However, if a large computational process is used to model the 
complex welding process, the thermophysical material properties 


on which the analysis is based must be accurately known. 


The axisymmetric finite element program which was developed 
by Muraki for analysis of thermal stress and metal movenent 
during welding is promising in that it showed good qualitative 
agreement with carefully generated experimental data. Heat 
flow and resulting temperature distributions within weldments 
can be accurately predicted by existing methods. Qualitative 
and quantitative comparison between analytical and measured 


temperature distributions showed remarkable coincidence. 
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СНАРТЕК УТТТ 


RECOMMENDATIONS 


Inasmuch as analytical methods to predict stresses and 
distortion resulting from thermal cycles can only be as good 
as the thermophysical material properties on which they are 
Based, it is recommended that studies be undertaken to 
accurately determine the temperature dependence of high 
strength steel properties at alevated temperatures. This 
would enhance further efforts to simulate the welding 


process. 


For further study of computer simulation of multi-pass 
welding, it is recommended that a micro-processor be used 
to process welding parameter data. Reduction of strains 
requiring temperature compensation is a tedious task, 
especially if measurements are made at several locations 
and recorded on one photo tape. This additional source of 
error could be reduced or eliminated if parameters of 
interest could be calculated internally and displayed 


и ресъту. 


High temperature strain gages could significantly con- 
tribute toward a more meaningful study. The 400 °F limitation 
of gages used in this study severely restricted their loca- 


Eion, 
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